Genomic characters of Anaplasma bovis and genetic diversity in China.

The emergence of Anaplasma bovis or A. bovis-like infection in humans from China and the United States of America has raised concern about the public health importance of this pathogen. Although A. bovis has been detected in a wide range of ticks and mammals in the world, no genome of the pathogen is available up to now, which has prohibited us from better understanding the genetic basis for its pathogenicity. Here we describe an A. bovis genome from metagenomic sequencing of an infected goat in China. Anaplasma bovis had the smallest genome of the genus Anaplasma, and relatively lower GC content. Phylogenetic analysis of single-copy orthologue sequence showed that A. bovis was closely related to A. platys and A. phagocytophilum, but relatively far from intraerythrocytic Anaplasma species. Anaplasma bovis had 116 unique orthogroups and lacked 51 orthogroups in comparison to other Anaplasma species. The virulence factors of A. bovis were significantly less than those of A. phagocytophilum, suggesting less pathogenicity of A. bovis. When tested by specific PCR assays, A. bovis was detected in 23 of 29 goats, with an infection rate up to 79.3% (95% CI: 64.6% ∼94.1%). The phylogenetic analyses based on partial 16S rRNA, gltA and groEL genes indicated that A. bovis had high genetic diversity. The findings of this study lay a foundation for further understanding of the biological characteristics and genetic diversity of A. bovis, and will facilitate the formulation of prevention and control strategies.

Characterization of a pangolin SARS-CoV-2-related virus isolate that uses the human ACE2 receptor.

Various SARS-CoV-2-related coronaviruses have been increasingly identified in pangolins, showing a potential threat to humans. Here we report the infectivity and pathogenicity of the SARS-CoV-2-related virus, PCoV-GX/P2V, which was isolated from a Malayan pangolin (Manis javanica). PCoV-GX/P2V could grow in human hepatoma, colorectal adenocarcinoma cells, and human primary nasal epithelial cells. It replicated more efficiently in cells expressing human angiotensin-converting enzyme 2 (hACE2) as SARS-CoV-2 did. After intranasal inoculation to the hACE2-transgenic mice, PCoV-GX/P2V not only replicated in nasal turbinate and lungs, but also caused interstitial pneumonia, characterized by infiltration of mixed inflammatory cells and multifocal alveolar hemorrhage. Existing population immunity established by SARS-CoV-2 infection and vaccination may not protect people from PCoV-GX/P2V infection. These findings further verify the hACE2 utility of PCoV-GX/P2V by in vivo experiments using authentic viruses and highlight the importance for intensive surveillance to prevent possible cross-species transmission.

An integrated data analysis reveals distribution, hosts, and pathogen diversity of Haemaphysalis concinna.

BACKGROUND: Haemaphysalis concinna, carrying multiple pathogens, has attracted increasing attention because of its expanded geographical range and significant role in disease transmission. This study aimed to identify the potential public health risks posed by H. concinna and H. concinna-associated pathogens. METHODS: A comprehensive database integrating a field survey, literature review, reference book, and relevant websites was developed. The geographical distribution of H. concinna and its associated pathogens was illustrated using ArcGIS. Meta-analysis was performed to estimate the prevalence of H. concinna-associated microbes. Phylogenetic and geographical methods were used to investigate the role of birds in the transmission of H. concinna-associated microbes. The potential global distribution of H. concinna was predicted by ecological niche modeling. RESULTS: Haemaphysalis concinna was distributed in 34 countries across the Eurasian continent, predominantly in China, Russia, and Central Europe. The tick species carried at least 40 human pathogens, including six species in the Anaplasmataceae family, five species of Babesia, four genospecies in the complex Borrelia burgdorferi sensu lato, ten species of spotted fever group rickettsiae, ten species of viruses, as well as Francisella, Coxiella, and other bacteria. Haemaphysalis concinna could parasitize 119 host species, with nearly half of them being birds, which played a crucial role in the long-distance transmission of tick-borne microbes. Our predictive modeling suggested that H. concinna could potentially survive in regions where the tick has never been previously recorded such as central North America, southern South America, southeast Oceania, and southern Africa. CONCLUSIONS: Our study revealed the wide distribution, broad host range, and pathogen diversity of H. concinna. Authorities, healthcare professionals, and the entire community should address the growing threat of H. concinna and associated pathogens. Tick monitoring and control, pathogen identification, diagnostic tools, and continuous research should be enhanced.

The complete mitogenome of Argas vulgaris (Filippova, 1961) and its phylogenetic status in subgenus Argas (Acari: Argasidae).

Soft ticks (Ixodida: Argasidae) are ectoparasites of terrestrial vertebrates with worldwide distributions. As one representative group of Argasidae, the genus Argas has an important vectorial role in transmitting zoonotic diseases. However, our knowledge of the subgenus Argas in China is still limited, as most literature only lists occurrence records or describes specific case reports without providing detailed morphological characteristics and further molecular data. This study aims to characterize Argas vulgaris through complete mitochondrial sequencing and morphological diagnostic techniques based on a batch of adult specimens collected from Ningxia Hui Autonomous Regions (NXHAR), North China. The morphology and microstructures of Ar. vulgaris and other lectotypes of argasid ticks in the subgenus Argas were also observed using a stereomicroscope. Following DNA extraction and sequencing, a complete mitochondrial sequence of Ar. vulgaris was assembled and analyzed within a phylogenetic context. The 14,479 bp mitogenome of Ar. vulgaris consists of 37 genes, including 13 genes for protein coding, two for ribosomal RNA, 22 for transfer RNA, and one for control region (D-loops). Phylogenetic analysis of Ar. vulgaris showed 98.27%-100% nucleotide identity with Ar. japonicus, indicating a close relationship between the two tick species. The morphological diagnostic features to differentiate Ar. vulgaris from other ticks within the subgenus Argas included the location of the anus and setae on the anterior lip of the female genital aperture. This study provided high-resolution scanning electron microscope images of female Ar. vulgaris and corresponding molecular data, representing valuable resources for future accurate species identification.

Isolation and characterization of a pangolin-borne HKU4-related coronavirus that potentially infects human-DPP4-transgenic mice.

We recently detected a HKU4-related coronavirus in subgenus Merbecovirus (named pangolin-CoV-HKU4-P251T) from a Malayan pangolin1. Here we report isolation and characterization of pangolin-CoV-HKU4-P251T, the genome sequence of which is closest to that of a coronavirus from the greater bamboo bat (Tylonycteris robustula) in Yunnan Province, China, with a 94.3% nucleotide identity. Pangolin-CoV-HKU4-P251T is able to infect human cell lines, and replicates more efficiently in cells that express human-dipeptidyl-peptidase-4 (hDPP4)-expressing and pangolin-DPP4-expressing cells than in bat-DPP4-expressing cells. After intranasal inoculation with pangolin-CoV-HKU4-P251, hDPP4-transgenic female mice are likely infected, showing persistent viral RNA copy numbers in the lungs. Progressive interstitial pneumonia developed in the infected mice, characterized by the accumulation of macrophages, and increase of antiviral cytokines, proinflammatory cytokines, and chemokines in lung tissues. These findings suggest that the pangolin-borne HKU4-related coronavirus has a potential for emerging as a human pathogen by using hDPP4.

Virome diversity shaped by genetic evolution and ecological landscape of Haemaphysalis longicornis.

BACKGROUND: Haemaphysalis longicornis is drawing attentions for its geographic invasion, extending population, and emerging disease threat. However, there are still substantial gaps in our knowledge of viral composition in relation to genetic diversity of H. longicornis and ecological factors, which are important for us to understand interactions between virus and vector, as well as between vector and ecological elements. RESULTS: We conducted the meta-transcriptomic sequencing of 136 pools of H. longicornis and identified 508 RNA viruses of 48 viral species, 22 of which have never been reported. Phylogenetic analysis of mitochondrion sequences divided the ticks into two genetic clades, each of which was geographically clustered and significantly associated with ecological factors, including altitude, precipitation, and normalized difference vegetation index. The two clades showed significant difference in virome diversity and shared about one fifth number of viral species that might have evolved to "generalists." Notably, Bandavirus dabieense, the pathogen of severe fever with thrombocytopenia syndrome was only detected in ticks of clade 1, and half number of clade 2-specific viruses were aquatic-animal-associated. CONCLUSIONS: These findings highlight that the virome diversity is shaped by internal genetic evolution and external ecological landscape of H. longicornis and provide the new foundation for promoting the studies on virus-vector-ecology interaction and eventually for evaluating the risk of H. longicornis for transmitting the viruses to humans and animals. Video Abstract.

Investigating the pathogens associated with Dermacentor nuttalli and its global distribution: A study integrating metagenomic sequencing, meta-analysis and niche modeling.

Dermacentor nuttalli, a member of family Ixodidae and genus Dermacentor, is predominantly found in North Asia. It transmits various pathogens of human and animal diseases, such as Lymphocytic choriomeningitis mammarenavirus and Brucella ovis, leading to severe symptoms in patients and posing serious hazards to livestock husbandry. To profile pathogen abundances of wild D. nuttalli, metagenomic sequencing was performed of four field-collected tick samples, revealing that Rickettsia, Streptomyces, and Pseudomonas were the most abundant bacterial genera in D. nuttalli. Specifically, four nearly complete Rickettsia genomes were assembled, closely relative to Rickettsia conorii subsp. raoultii. Then, a comprehensive meta-analysis was performed to evaluate its potential threats based on detected pathogens and geographical distribution positions reported in literature, reference books, related websites, and field surveys. At least 48 pathogens were identified, including 20 species of bacteria, seven species of eukaryota, and 21 species of virus. Notably, Rickettsia conorii subsp. raoultii, Coxiella burnetii, and Brucella ovis displayed remarkably high positivity rates, which were known to cause infectious diseases in both humans and livestock. Currently, the primary distribution of D. nuttalli spans China, Mongolia, and Russia. However, an additional 14 countries in Asia and America that may also be affected by D. nuttalli were identified in our niche model, despite no previous reports of its presence in these areas. This study provides comprehensive data and analysis on the pathogens carried by D. nuttalli, along with documented and potential distribution, suggesting an emerging threat to public health and animal husbandry. Therefore, there is a need for heightened surveillance and thorough investigation of D. nuttalli.

Changes in Cytokine Levels in Patients with Severe Fever with Thrombocytopenia Syndrome Virus.

Purpose: To characterize the cytokine profile of patients with severe fever with thrombocytopenia syndrome (SFTS) in relation to disease severity. Patients and Methods: 60 laboratory-confirmed SFTS patients and 12 healthy individuals from multi-centers in Shandong Province of China were included, and all patients were divided into fatal patients (9) and recovered patients (51) due to their final outcomes. Multiplex-microbead immunoassays were conducted to estimate levels of 27 cytokines in the sera of patients and controls. Results: The results showed that levels of IL-2, IL-4, IL-6, IL-7, IL-8, IL-15, IL-1RA, G-CSF, GM-CSF, IFN-γ, TNF-α, basic FGF, PDGF-BB, RANTES, IP-10, MIP-1α, MIP-1ß, MCP-1, and Eotaxin differed significantly among the SFTS fatal patients, recovered patients, and the healthy controls (all p<0.05). Compared to the healthy controls, the fatal patients and recovered patients had reduced levels of IL-2, IL-4, IL-7, PDGF-BB, RANTES, and Eotaxin, while the levels of PDGF-BB and RANTES were significantly lower in fatal patients compared to recovered patients. The increasing levels of IL-6, IL-8, IL-15, IL-1RA, G-CSF, GM-CSF, IFN-γ, TNF-α, basic FGF, IP-10, MIP-1α, MIP-1ß, and MCP-1 were observed in fatal patients (all p<0.05), and the levels of IL-6, IP-10, MIP-1α, and MCP-1 were significantly higher than other two groups. The Spearman correlation analysis indicated a positive correlation between platelet count and PDGF-BB levels (p<0.05), while the white blood cell count had a negative correlation with MIP-1 level (p<0.05). Conclusion: The research exhibited that the SFTS virus (SFTSV) caused an atypical manifestation of cytokines. The levels of IL-6, IP-10, MIP-1α, and MCP-1 had been observed a positive association with the severity of the illness.

Comorbidity increases the risk of pulmonary tuberculosis: a nested case-control study using multi-source big data.

BACKGROUND: Some medical conditions may increase the risk of developing pulmonary tuberculosis (PTB); however, no systematic study on PTB-associated comorbidities and comorbidity clusters has been undertaken. METHODS: A nested case-control study was conducted from 2013 to 2017 using multi-source big data. We defined cases as patients with incident PTB, and we matched each case with four event-free controls using propensity score matching (PSM). Comorbidities diagnosed prior to PTB were defined with the International Classification of Diseases-10 (ICD-10). The longitudinal relationships between multimorbidity burden and PTB were analyzed using a generalized estimating equation. The associations between PTB and 30 comorbidities were examined using conditional logistic regression, and the comorbidity clusters were identified using network analysis. RESULTS: A total of 4265 cases and 17,060 controls were enrolled during the study period. A total of 849 (19.91%) cases and 1141 (6.69%) controls were multimorbid before the index date. Having 1, 2, and ≥ 3 comorbidities was associated with an increased risk of PTB (aOR 2.85-5.16). Fourteen out of thirty comorbidities were significantly associated with PTB (aOR 1.28-7.27), and the associations differed by sex and age. Network analysis identified three major clusters, mainly in the respiratory, circulatory, and endocrine/metabolic systems, in PTB cases. CONCLUSIONS: Certain comorbidities involving multiple systems may significantly increase the risk of PTB. Enhanced awareness and surveillance of comorbidity are warranted to ensure early prevention and timely control of PTB.

Epidemiological behaviour and interventions of malaria in Niger, 2010-2019: a time-series analysis of national surveillance data.

BACKGROUND: Malaria remains a significant public health concern in Niger, with the number of cases increasing from 592,334 in 2000 to 3,138,696 in 2010. In response, a concerted campaign against the disease has been initiated. However, the implementation of these malaria interventions and their association with epidemiological behaviour remains unclear. METHODS: A time-series study was conducted in Niger from 2010 to 2019. Multiple data sources concerning malaria were integrated, encompassing national surveillance data, Statistic Yearbook, targeted malaria control interventions, and meteorological data. Incidence rate, mortality rate, and case fatality ratio (CFR) by different regions and age groups were analysed. Joinpoint regression models were used to estimate annual changes in malaria. The changes in coverage of malaria interventions were evaluated. RESULTS: Between 2010 to 2019, the incidence rate of malaria decreased from 249.43 to 187.00 cases per 1,000 population in Niger. Niamey had a high annual mean incidence rate and the lowest CFR, while Agadez was on the contrary. Joinpoint regression analysis revealed a declining trend in malaria incidence for all age groups except the 10-24 years group, and the mortality rate and the CFR initially decreased followed by an increase in all age groups. Niger has implemented a series of malaria interventions, with the major ones being scaled up to larger populations during the study period. CONCLUSIONS: The scale-up of multi-interventions in Niger has significantly reduced malaria incidence, but the rise in mortality rate and CFR addresses the challenges in malaria control and elimination. Malaria endemic countries should enhance surveillance of malaria cases and drug resistance in Plasmodium, improve diagnosis and treatment, expand the population coverage of insecticide-treated bed nets and seasonal malaria chemoprevention, and strengthen the management of severe malaria cases.

Systematic review and integrated data analysis reveal diverse pangolin-associated microbes with infection potential.

There has been increasing global concern about the spillover transmission of pangolin-associated microbes. To assess the risk of these microbes for emergence as human pathogens, we integrated data from multiple sources to describe the distribution and spectrum of microbes harbored by pangolins. Wild and trafficked pangolins have been mainly recorded in Asia and Africa, while captive pangolins have been reported in European and North American countries. A total of 128 microbes, including 92 viruses, 25 bacteria, eight protists, and three uncharacterized microbes, have been identified in five pangolin species. Out of 128 pangolin-associated microbes, 31 (including 13 viruses, 15 bacteria, and three protists) have been reported in humans, and 54 are animal-associated viruses. The phylogenetic analysis of human-associated viruses carried by pangolins reveals that they are genetically close to those naturally circulating among human populations in the world. Pangolins harbor diverse microbes, many of which have been previously reported in humans and animals. Abundant viruses initially detected in pangolins might exhibit risks for spillover transmission.

Epidemiological and phylogenetic characteristics of emerging Anaplasma capra: A systematic review with modeling analysis.

Anaplasma capra, an emerging tick-borne pathogen, has caused a lot of concern since initially recognized in goats and patients in China in 2015, and has been reported in a wide range of domestic and wild animals as well as ticks worldwide, posing a threat to public health. In this systematic review, we established a comprehensive database to acquire the distribution and prevalence status of this pathogen, and collected all sequences of A. capra to summarize the details of genetic diversity by phylogenetic analysis. We also predicted the possible global distribution of A. capra by using ecological niche model. A. capra has been known to distribute in 18 countries across Asia, Europe and Africa. A total of 19 species of mammals from seven families have been reported as hosts, and domestic ruminants including goats, sheep and cattle were the major hosts. At least 8 tick species of 4 genera have been reported to carry A. capra, and Haemaphysalis longicornis was most commonly infected. Sheep and Rhipicephalus microplus had the highest positive rates among animals and ticks. Phylogenetic analysis based on gltA and groEL genes revealed that A. capra could primarily be divided into two clusters related to geographic location and animal hosts. The predictive model showed that the most suitable habitats for presence of A. capra were mainly located in Asia and eastern Europe. These cumulative data regarding A. capra of our study lay a foundation for the subsequent exploration of this emerging tick-borne pathogen.

Genomic Characteristics of Emerging Intraerythrocytic Anaplasma capra and High Prevalence in Goats, China.

Anaplasma capra is an emerging tickborne human pathogen initially recognized in China in 2015; it has been reported in ticks and in a wide range of domestic and wild animals worldwide. We describe whole-genome sequences of 2 A. capra strains from metagenomic sequencing of purified erythrocytes from infected goats in China. The genome of A. capra was the smallest among members of the genus Anaplasma. The genomes of the 2 A. capra strains contained comparable G+C content and numbers of pseudogenes with intraerythrocytic Anaplasma species. The 2 A. capra strains had 54 unique genes. The prevalence of A. capra was high among goats in the 2 endemic areas. Phylogenetic analyses revealed that the A. capra strains detected in this study were basically classified into 2 subclusters with those previously detected in Asia. Our findings clarify details of the genomic characteristics of A. capra and shed light on its genetic diversity.

Successful treatment with doxycycline monotherapy for human infection with Babesia venatorum (Babesiidae, Sporozoa) in China: a case report and proposal for a clinical regimen.

BACKGROUND: Human babesiosis is a worldwide disease caused by intraerythrocytic protozoa of the genus Babesia. It is transmitted by bites from ixodid ticks, and mechanically transmitted by blood transfusion. It is primarily treated with quinine and/or atovaquone, which are not readily available in China. In this study, we developed a novel treatment regimen involving doxycycline monotherapy in a patient with severe Babesia venatorum infection as an alternative therapeutic medication. The aim of our study is to provide a guidance for clinical practice treatment of human babesiosis. CASE PRESENTATION: A 73-year-old man who had undergone splenectomy and blood transfusion 8 years prior, presented with an unexplained fever, headache, and thrombocytopenia, and was admitted to the Fifth Medical Center of the PLA General Hospital. He was diagnosed with B. venatorum infection by morphological review of thin peripheral blood smears, which was confirmed by multi-gene polymerase chain reaction (PCR), and sequencing of the entire 18s rRNA and partial ß-tubulin encoding genes, as well as isolation by animal inoculation. The doxycycline monotherapy regimen (peros, 0.1 g bisindie) was administered following pharmacological guidance and an effective outcome was observed. The patient recovered rapidly following the doxycycline monotherapy. The protozoan load in peripheral blood samples decreased by 88% in hematocrit counts after 8 days, and negative PCR results were obtained after 90 days of follow-up at the hospital. The treatment lasted for 3 months without any side effects or sequelae. The nine-month follow-up survey of the patient did not reveal any signs of recrudescence or anti-babesial tolerance. CONCLUSIONS: We have reported a clinical case of successful doxycycline monotherapy for human babesiosis caused by B. venatorum, which provides an optional medical intervention for human babesiosis.

Genomic Characterization of Theileria luwenshuni Strain Cheeloo.

Theileria, a tick-borne intracellular protozoan, can cause infections of various livestock and wildlife around the world, posing a threat to veterinary health. Although more and more Theileria species have been identified, genomes have been available only from four Theileria species to date. Here, we assembled a whole genome of Theileria luwenshuni, an emerging Theileria, through next-generation sequencing of purified erythrocytes from the blood of a naturally infected goat. We designated it T. luwenshuni str. Cheeloo because its genome was assembled by the researchers at Cheeloo College of Medicine, Shandong University, China. The genome of T. lunwenshuni str. Cheeloo was the smallest in comparison with the other four Theileria species. T. luwenshuni str. Cheeloo possessed the fewest gene gains and gene family expansion. The protein count of each category was always comparable between T. luwenshuni str. Cheeloo and T. orientalis str. Shintoku in the Eukaryote Orthologs annotation, though there were remarkable differences in genome size. T. luwenshuni str. Cheeloo had lower counts than the other four Theileria species in most categories at level 3 of Gene Ontology annotation. Kyoto Encyclopedia of Genes and Genomes annotation revealed a loss of the c-Myb in T. luwenshuni str. Cheeloo. The infection rate of T. luwenshuni str. Cheeloo was up to 81.5% in a total of 54 goats from three flocks. The phylogenetic analyses based on both 18S rRNA and cox1 genes indicated that T. luwenshuni had relatively low diversity. The first characterization of the T. luwenshuni genome will promote better understanding of the emerging Theileria. IMPORTANCE Theileria has led to substantial economic losses in animal husbandry. Whole-genome sequencing data of the genus Theileria are currently limited, which has prohibited us from further understanding their molecular features. This work depicted whole-genome sequences of T. luwenshuni str. Cheeloo, an emerging Theileria species, and reported a high prevalence of T. luwenshuni str. Cheeloo infection in goats. The first assembly and characterization of T. luwenshuni genome will benefit exploring the infective and pathogenic mechanisms of the emerging Theileria to provide scientific basis for future control strategies of theileriosis.

Continuous evolution and emerging lineage of seasonal human coronaviruses: A multicenter surveillance study.

The seasonal human coronaviruses (HCoVs) have zoonotic origins, repeated infections, and global transmission. The objectives of this study are to elaborate the epidemiological and evolutionary characteristics of HCoVs from patients with acute respiratory illness. We conducted a multicenter surveillance at 36 sentinel hospitals of Beijing Metropolis, China, during 2016-2019. Patients with influenza-like illness (ILI) and severe acute respiratory infection (SARI) were included, and submitted respiratory samples for screening HCoVs by multiplex real-time reverse transcription-polymerase chain reaction assays. All the positive samples were used for metatranscriptomic sequencing to get whole genomes of HCoVs for genetical and evolutionary analyses. Totally, 321 of 15 677 patients with ILI or SARI were found to be positive for HCoVs, with an infection rate of 2.0% (95% confidence interval, 1.8%-2.3%). HCoV-229E, HCoV-NL63, HCoV-OC43, and HCoV-HKU1 infections accounted for 18.7%, 38.3%, 40.5%, and 2.5%, respectively. In comparison to ILI cases, SARI cases were significantly older, more likely caused by HCoV-229E and HCoV-OC43, and more often co-infected with other respiratory pathogens. A total of 179 full genome sequences of HCoVs were obtained from 321 positive patients. The phylogenetical analyses revealed that HCoV-229E, HCoV-NL63 and HCoV-OC43 continuously yielded novel lineages, respectively. The nonsynonymous to synonymous ratio of all key genes in each HCoV was less than one, indicating that all four HCoVs were under negative selection pressure. Multiple substitution modes were observed in spike glycoprotein among the four HCoVs. Our findings highlight the importance of enhancing surveillance on HCoVs, and imply that more variants might occur in the future.

Geographical distribution of Ixodes persulcatus and associated pathogens: Analysis of integrated data from a China field survey and global published data.

The increasing incidence and range expansion of tick-borne diseases have caused global threats to human and animal health under the background of climate and socioeconomic changes. As an efficient vector in transmission of tick-borne diseases, a growing burden caused by Ixodes persulcatus and associated pathogens could not be underestimated. This study summarized the distribution, hosts, and pathogens of I. persulcatus, and predicted the suitable habitats of this tick species worldwide. An integrated database involving a field survey, reference book, literature review, and related website was constructed. Location records of I. persulcatus and associated pathogens were incorporated into distribution maps using ArcGIS software. Positive rates for I. persulcatus-associated agents were estimated by meta-analysis. The global distribution of the tick species was predicted using Maxent model. I. persulcatus was distributed in 14 countries across the Eurasian continent, involving Russia, China, Japan, and several Baltic Sea states, which ranged between 21°N to 66°N. The tick species fed on 46 species of hosts, and 51 tick-borne agents could be harbored by I. persulcatus. The predictive model showed that I. persulcatus could be predominantly distributed in northern Europe, western Russia, and northern China. Our study fully clarified the potential public health risks posed by I. persulcatus and I. persulcatus-borne pathogens. Surveillance and control measures of tick-borne diseases should be enhanced to promote the health of humans, animals, and ecosystems.

New insights into the impact of microbiome on horizontal and vertical transmission of a tick-borne pathogen.

BACKGROUND: The impact of host skin microbiome on horizontal transmission of tick-borne pathogens , and of pathogen associated transstadial and transovarial changes in tick microbiome are largely unknown, but are important to control increasingly emerging tick-borne diseases worldwide. METHODS: Focusing on a rickettsiosis pathogen, Rickettsia raoultii, we used R. raoultii-positive and R. raoultii-negative Dermacentor spp. tick colonies to study the involvement of skin microbiota in cutaneous infection with rickettsiae in laboratory mice, and the function of the tick microbiome on maintenance of rickettsiae through all tick developmental stages (eggs, larvae, nymphs, adults) over two generations. RESULTS: We observed changes in the skin bacteria community, such as Chlamydia, not only associated with rickettsial colonization but also with tick feeding on skin. The diversity of skin microbiome differed between paired tick-bitten and un-bitten sites. For vertical transmission, significant differences in the tick microbiota between pathogenic rickettsia-positive and -negative tick chorts was observed across all developmental stages at least over two generations, which appeared to be a common pattern not only for R. raoultii but also for another pathogenic species, Candidatus Rickettsia tarasevichiae. More importantly, bacterial differences were complemented by functional shifts primed for genetic information processing during blood feeding. Specifically, the differences in tick microbiome gene repertoire between pathogenic Rickettsia-positive and -negative progenies were enriched in pathways associated with metabolism and hormone signals during vertical transmission. CONCLUSIONS: We demonstrate that host skin microbiome might be a new factor determining the transmission of rickettsial pathogens through ticks. While pathogenic rickettsiae infect vertebrate hosts during blood-feeding by the tick, they may also manipulate the maturation of the tick through changing the functional potential of its microbiota over the tick's life stages. The findings here might spur the development of new-generation control methods for ticks and tick-borne pathogens. Video Abstract.

Metavirome of 31 tick species provides a compendium of 1,801 RNA virus genomes.

The increasing prevalence and expanding distribution of tick-borne viruses globally have raised health concerns, but the full repertoire of the tick virome has not been assessed. We sequenced the meta-transcriptomes of 31 different tick species in the Ixodidae and Argasidae families from across mainland China, and identified 724 RNA viruses with distinctive virome compositions among genera. A total of 1,801 assembled and complete or nearly complete viral genomes revealed an extensive diversity of genome architectures of tick-associated viruses, highlighting ticks as a reservoir of RNA viruses. We examined the phylogenies of different virus families to investigate virome evolution and found that the most diverse tick-associated viruses are positive-strand RNA virus families that demonstrate more ancient divergence than other arboviruses. Tick-specific viruses are often associated with only a few tick species, whereas virus clades that can infect vertebrates are found in a wider range of tick species. We hypothesize that tick viruses can exhibit both 'specialist' and 'generalist' evolutionary trends. We hope that our virome dataset will enable much-needed research on vertebrate-pathogenic tick-associated viruses.